Spin gap and Luttinger liquid description of the NMR relaxation in carbon nanotubes 
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Recent NMR experiments by Singer et al. [Singer et al. Phys. Rev. Lett. 95, 236403 (2005).] 
showed a deviation from Fermi-liquid behavior in carbon nanotubes with an energy gap evident at 
low temperatures. Here, a comprehensive theory for the magnetic field and temperature dependent 
NMR spin-lattice relaxation is given in the framework of the Tomonaga-Luttinger liquid. The 
low temperature properties are governed by a gapped relaxation due to a spin gap (~ 30 K), which 
crosses over smoothly to the Luttinger liquid behaviour with increasing temperature. 
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PACS numbers: 71.10.Pm,71.20.Tx,76.60.-k 



Low dimensional carbonaceous systems, fuUerenes, 
carbon nanotubes (CNTs), and graphene display a 
rich variety of exotic states and strongly correlated 
phenomena. These include superconductivity in al- 
kali doped fuUerenesi, quantized transport in single- 
wall carbon nanotubes (SWCNTs), and massless Dirac 
quasi-particles showing a halve integer quantum Hall- 
effect in graphene even at room temperature. A com- 
pelling correlated state of one-dimensional systems is the 
Tomonaga-Luttinger liquid (TLL) state. The TLL state 
has been suggested to describe the low energy proper- 
ties of CNTs with a single shell, the single-wall car- 
bon nanotubes^ii^i^i^i^. Transport^ and photoemission 
studies^'^ provided evidence for the existence of the TLL 
state in SWCNTs. In these studies, power-law behavior 
of temperature and bias dependent conductivity and a 
power-law Fermi edge was observed, respectively. 

Nuclear magnetic resonance (NMR) is a powerful 
method to characterize correlated states of materials as 
it is sensitive to the density of states near the Fermi 
edge. For a material with a Fermi-liquid state, the tem- 
perature dependent spin-lattice relaxation time, Ti fol- 
lows the so-called Korringa temperature dependence for 
which 1/TiT is constant. Recently, ^^C enriched SWC- 
NTs were grown inside carbon nanotubes from ^^C en- 
riched fullerenesiS. This allowed a high precision mea- 
surement of Ti in small diameter SWCNTs by Singer et 
aliii. A tentative fit of the experiments with a gapped 
Fermi liquid type density of states (DOS) indicated over- 
all agreement but obvious discrepancies in detail. When 
the magnetic field and temperature dependent data for 
Ti were fitted with this phenomenological model a gap 
at the Fermi surface with 2A ~ 40 K opened already 
above room temperature, i.e. its Tc is larger than 300 K. 
This strongly violates the 2A/kBT > 3.52 relation, thus 
simple mean field theories are not applicable--. Also, 
the phenomenological description can not account for the 
strong overshoot of 1/TiT when T approaches the gap. 

Here, we analyze the NMR results in the framework 
of the Luttinger liquid and Luther-Emery liquid pictures 



(interacting one-dimensional electrons without and with 
a gap, respectively). At high temperatures, the former 
dominates, while the latter accounts for the dominance 
of a spin gap at low temperatures. We show that the 
temperature and magnetic field dependent Ti of ^^C can 
be explained in the TLL scenario with a small spin gap 
of the order of 30 K. 

The most striking difference between Fermi and Lut- 
tinger liquids is the anomalous power-law dependence of 
correlation functions at low energies in the latter (includ- 
ing the NMR relaxation rate), with critical exponents 
changing continuously with the interaction strength. 
This is to be contrasted to a Fermi liquid, where the 
quasi particle picture (electron) holds, implying critical 
exponents fixed to an integer. In a Luttinger liquid, low 
energy properties are described by critical phenomena of 
collective modes. 

In the following, we construct the low energy theory 
for CNT with the use of bosonization. We determine 
the various contributions to the NMR relaxation rate at 
high and low temperatures, and fit the obtained formula 
to the experimental results. 

The Hamiltonian describing the kinetic energy of elec- 
trons on SWNT is given by^iii^i^ 



dxY,r^tc.,{x)d,^r<.„{x), (1) 



where only two transport bands (a = ±) are taken into 
account, and for each band, there is a right and left mov- 
ing branch (r = ±) with linear dispersion, cr stands for 
the spin. The forward scattering interaction reads as 



H 



FS 



dxdyp{x)V{x - y)p{y), 



(2) 



where p{x) = J2raa'^tacrix)'^raaix), V (x) IS the 

Coulomb potential. For the moment, we neglect the umk- 
lapp and backward scattering processes. The Hamilto- 
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nian is bosonized via 



exp [i{qFr + kFa)x + i{r6aa + 4>aa)] 



(3) 

where kp denotes the undoped Fermi surface of a 
graphene layer, qp determines the band fiUingS and 
^ kp to ensure the apphcabihty of the low en- 
ergy continuum approximation, rjrao^s are Majorana 
fermions, standing for the proper anticommutation re- 
lations between different branches. The phase fields are 
expressed by the symmetric and antisymmetric charge 
and spin modes [5 — ±) between the valleys of the charge 
and spin excitations (j — c, s) as 
aOc- + aaOs- and 2(j)aa = 4>c+ + o" 
This leads to the Hamiltonian 



(tOs+ + 



Ho — Hkin 



H 



FS 



dx 



(4) 



Vjs — vp/Kjs- For long range Coulomb interaction, 
Kc+ < 1 and is usually around 0.2, while for the 
other modes, Kjs > 1, but only slightly greater than 

Other interaction terms such as umklapp and backscat- 
tering have exhaustively been studied in Refs. and 
we refrain from their analyzis here. These parallel closely 
to the investigation of coupled Luttinger liquid chains. 
Their effect can be summarized as follows: at half filling, 
umklapp scattering generates a finite charge (Mott) gap 
in both the symmetric and antisymmetric charge sector, 
while spin backscattering leads to the opening of a spin 
gap in the spin sectorsMi^. This is analogous to the spin 
ladder magnon excitations^'* found in the strong antifer- 
romagnetic coupling limit. Thus, the ground state at half 
filling is a Mott insulator with all excitations gappedi^. 
Away from half filling, the c-l- sector is gapless as dis- 
cussed before, while all other sectors remain massive. 
This means, that at low temperatures (T ^ A), expo- 
nentially activated behaviour is expected'^'^, while in the 
high temperature regime (T ^ A), the usual Luttinger 
liquid behaviour is restored. Since the strength of back- 
ward and umklapp interactions is inversely proportional 
to the diameter of the tube, these become important (and 
generate a gap) o nly for narrow tubes^i^ii^, such as the 
one in Refs. llOlflll with an inner tube diameter of 0.7 nm. 

We start our analysis in the high temperature region 
to characterize the TLL behaviour first. Hence, we post- 
pone the discussion of the gaps for the moment, and con- 
centrate on the temperature exponents in the NMR re- 
laxation rate, l/Ti. In general, 1/Ti measures the local 
dynamics of the spins, and it is related to the transverse 
spin susceptibility, which in the case of nanotubes, reads 
as 



X±(*^n) = - j dTe'^-^TrS+{x - G,T)5-(a; = 0,0)), 



(5) 



where S^{x) = Y.rar'a'^ta±i.x)^r'a'^{x). From this, 
after analytic continuation to real frequencies, the NMR 
relaxation time is 



(TiT)- 



lim 



(6) 



In general, there are four different length scales^ in 
S^{x), and the corresponding wavevectors are 1. qi = 0, 
2. q2 = 2kp, 3. (73 — 2qp and 4. 54 — 2{kp ± qp). These 
show up in the term to S^{x) at 1. a — a', r = r' , 2. 
a = —a', r = r', 3. a = a', r = —r' and 4. a = —a', 
r = — r', respectively. 
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exponent (7^) 
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a = a' , r — r' 
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a = —a', r — r' 
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a — a' , r — — r' 
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2{kF ± Qf) 


Q = ^ct , r = — r' 





TABLE I: The long time exponents of the different contribu- 
tions to the transverse spin susceptibility (~ t ''''^). These 
also determine the temperature dependence of the NMR re- 
laxation rate as (TiT)'^ ~ y7i/2-2^ 

In all cases, the corresponding contribution to the sus- 
ceptibility decays at long times as ~ l/f'/^ with the 
exponents listed in Table HI Based on Ref. 2, these 
can be approximated by 71 ~ 72 ~ 2 + 1 /Kg + Ks and 
73 « 74 « ifc + 2 + II Ks with Ks = Ks+ and K^ = A'c+. 
In the most general case with broken spin rotational sym- 
metry, all these exponents usually differ from each other. 
However, from realistic models^i^i^, one can conjecture 
their values as Kg- = Kc- = 1. In the presence of ex- 
ternal magnetic field, which couples to the s+ field and 
breaks the spin rotational invariance, marginal interac- 
tions change the exponents^^, hence K's can differ from 
their zero field value. This follows from the exact Bethe- 
Ansatz solution of the Hubbard model in magnetic field- 
as well. Based on all this, a simple power counting deter- 
mines the temperature exponents in the NMR relaxation 
rate as ~ TT./2-2. 

Having determined the high temperature Luttinger liq- 
uid behaviour, we now turn to the evaluation of the 
response function at low temperatures where the gaps 
are dominating. In general, at T = 0, all the correla- 
tion functions containing a gapped field acquire a factor 
^ exp(— AVt^ -I- w^a;^), which suppresses exponentially 
the long time-long distance power-law Luttinger liquid 
behaviour^ii^. This picture is further corroborated by 
the existing exact solutions of the attractive Hubbard 
model along the Luther- Emery line {K — 1/2), where 
strongly interacting gapped bosons are mapped onto non- 
interacting fermionsi^ii^, i.e. the non-interacting mas- 
sive Thirring model^S. In some cases, this mapping also 
allows for an explicit evaluation of certain correlation 
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functions, such as the finite momentum transverse spin 
susceptibihty^'''^°. 

At finite temperatures, this simple picture needs to be 
modified. From the exact solution of gapped systems 
such as the sine-Gordon modeli^ii^, it is well established 
that different type of excitations contribute to higher 
frequencies and temperatures^^. The simplest ones are 
the solitons and antisolitons (kinks), which, when attract 
each other, form breathers. Their contribution can sys- 
tematically be analyzed based on the " Form Factor Boot- 
strap Approach"—!^, which uses the integrability of the 
underlying model (in the case of the Hubbard model, for 
example) to determine the matrix elements of the appro- 
priate operators between the ground state and excited 
states. In this simple exponential term cannot 

account for their variety. 

The explicit determination of these matrix elements 
is certainly beyond the scope of the present investiga- 
tion. Instead, we utilize our knowledge about these exci- 
tations. At temperatures comparable to A, they start to 
contribute significantly to the response fimctions. This 
situation can be mimicked by retaining our simple expo- 
nential factor in the correlation functions, and adjusting 
the gap at each temperature to give reasonable contribu- 
tion. In other words, this means a A(r) gap function. 
In spirit, this approach is identical to the self-consistent 
harmonic approximation^, where one replaces the com- 
plicated interaction term (such as the cosine in the sine- 
Gordon model) by simpler quadratic ones, and optimizes 
its coefficient to minimize the free energy of the systent^^. 
Such an approach for the sine- Gordon model leads to a 
mass which monotonically decreases with temperature, 
while above certain T ^ A(r = 0), the potential gener- 
ating a gap for a field is completely washed due to the 
fluctuations of the field itself^^. Above these tempera- 
tures, the mass practically disappears. Such an approach 
gives very accurate results for the dependence of the gap 
on model parameters at zero temperatur o^'^i^^ . 

At finite temperatures, contributions to the real time 
transverse spin correlator read as 

|sinh(7rTt)|^*/^ ^ 4 

(7) 

which reduces to the standard Luttinger liquid correlator 
in the absence of the gap. Here A is the sum of the 
existing gaps. At half filling, all sectors are gappedi^, 
hence A is their sum. Awayfrom half filling, which is 
presumably the case in Ref. Uil, the spin sectors are fully 
gapped^., and the A is dominated by their contribution. 
From this, by Fourier transformation, we can evaluate 
the frequency dependent retarded correlation functions, 
from which we can determine the NMR relaxation rate 
as 



where 



F{T) = exp I 1 -, — . . ,1, — ^ — — X 



27tT 2 



^ - 



iA{T) K 

2ttt ^ y 



1 - 



iA(T) K 

2'kt ^ y 



(9) 



and K — 7i/2, A is a proportionality factor, which is de- 
termined by the contribution of the respective fermionic 
fields and cannot be obtained by the method used here, 
T{x) and ^(2;) are the gamma and digamma functions, 
respectively. At high temperatures, Eq. ^ exhibits 
a T^~^ behaviourS^, as is expected from Luttinger 
liquidai^ii^,. Based on the possible values of Luttinger liq- 
uid parameters, the contributions with 73 and 74 are the 
most divergent ones at low temperatures due to the pres- 
ence of the c+ mode, hence we will concentrate on them 
in the followings. They lead io K ^l + {Kc + l/Ks)/2. 




2aTT 

V 



K 



(8) 



FIG. 1: (Color online) The fit of Eq. ([8| to the experimental 
data (open circles) of Ref. [ij at _ff = 3.6 T (red dashed 
line) and 9.3 T (blue line). The inset shows the deduced Ks 
exponents as a function of the field for Kc = 0.2 (black circles) 
and 0.05 (red squares). 

In Fig. [Tl we fit our Eq. ([8]) to the experimen- 
tal data on the inner wall of metallic double-wall car- 
bon nanotubes as a function of temperature at magnetic 
fields H — 3.6 T and 9.3 T. To increase the strength 
of the NMR signal, selective enrichment of the inner 
shells using ^^C isotope was performed, hence mainly 
the properties of the inner tube has been probed by the 
magnetic resonance measurements. Here, we take A = 
Ao(l — T/Tc)^/'^, which is very close to that found in the 
self consistent solution of the soliton energy in the sine- 
Gordon model2^. With this, we can nicely account for 
the various complicated excitations of our model. Obvi- 
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ously, we expect a smooth crossover of A as the tempera- 
ture increases through to a smaU value {A{T)/T ^ 1) 
rather than a sharp drop to zero, but for computational 
purposes, we use the above sharp form. Our results are 
robust with respect to variations of the gap function, e.g. 
the quality of the fitting remains the same by changing 
the exponent, 1/4, within the range 0.1 - 0.7, since only 
a finite gap a low temperatures (T <^ Tc) and a tiny gap 
at high temperatures is required (T ^ Tc). For the fits, 
we use Tc « 13 K and Aq « 32 K. Theoretically25, the 
critical temperature, at which the cosine potential practi- 
cally disappears from the sine-Gordon model, is roughly 
determined as Tc « A(0)/e, in perfect agreement with 
our fitting. Above Tc, the gap does not make itself felt 
any more. For the Luttinger liquid parameters, we as- 
sume Kc = 0.2, in accordance with Refs. i and i, and 
deduce Ks = 1.07 for H = 3.6 T and = 0.87 for 
H — 9.3 T, which vary with the field but stay close 
to unity as expected^"^. For the sake of completeness, 
we give the numerical prefactors, which are found to be 
A = 1.4 X 10^ K/s, while upon reinserting original units, 
kB^an/hv — 10"'' 1/K, and v can be taken roughly as 
vp. The latter falls into the same order of magnitude, 
if we identify a, as a rough estimate, to be of the or- 
der of the honeycomb lattice constant 2.46 A), and 
vp ~ 10^ m/s, which again certifies our approach. Since 
the NMR relaxation rate has been measured at 290 K for 
different fieldsii, we can also determine the field depen- 
dence of the Ks exponent, which is shown in the inset of 
Fig. [l]for two different, field independent values of Kc. 
It stays close to 1 and decreases with field. The other 
exponents Kc± and Ks- can also vary with the field, but 



since no exact solutions are available for the Luttinger 
liquid exponents of carbon nanotubes unlike for the Hub- 
bard modeli^, we choose the simplest possible approach 
of assigning all the field dependence to the s+ field. It 
is evident from the figure, that the larger Kc, the larger 
Kg, hence a suitable increase of the former in magnetic 
field can result in an increment of the latter. Given the 
simplicity of our scheme, the obtained fits are excellent 
and the parameters are reasonable. Finally, we mention 
the possibility, that upon doping, the properties of the 
c+ mode are expected to change (e.g. Kc+) similarly to 
the change of the s+ mode in a magnetic field. 

In conclusion, we have analyzed a recent experimental 
dataiJ of the NMR relaxation rate of CNT. With the use 
of bosonization, we find that the low temperature region 
is well described by a spin gapped Luther-Emery liquid. 
Relaxation is dominated by exponentially activated be- 
haviour with a gap ~ 30 K. This crosses over to the usual 
Luttinger liquid picture at high temperatures compared 
to the gap, with characteristic non-integer power-law ex- 
ponents. Finally, we predict, that the observed spin gap 
can significantly be reduced and even completely sup- 
pressed through a quantum phase transition by experi- 
mentally accessible magnetic fields {^-^ 25 — 30 T) due to 
its small value, similarly^'^'^'' to Cu2(C5Hi2N2)2Cl4. 
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